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Abstract 

Prompt D meson and non-prompt J/y/ yields are studied as a function of the multiplicity of charged 
particles produced in inelastic proton-proton collisions at a centre-of-mass energy of ^/s = 7 TeV. 
The results are reported as a ratio between yields in a given multiplicity interval normalised to the 
multiplicity-integrated ones (relative yields). They are shown as a function of the multiplicity of 
charged particles normalised to the average value for inelastic collisions (relative charged-particle 
multiplicity). D**, D+ and D*+ mesons are measured in five pj intervals from 1 GeV/c to 20 GeV/c 
and for |y| < 0.5 via their hadronic decays. The D-meson relative yield is found to increase with 
increasing charged-particle multiplicity. For events with multiplicity six times higher than the av¬ 
erage multiplicity of inelastic collisions, a yield enhancement of a factor about 15 relative to the 
multiplicity-integrated yield in inelastic collisions is observed. The yield enhancement is indepen¬ 
dent of transverse momentum within the uncertainties of the measurement. The D^'-meson relative 
yield is also measured as a function of the relative multiplicity at forward pseudo-rapidity. The 
non-prompt J/1//, i.e. the B hadron, contribution to the inclusive J/y/ production is measured in the 
di-electron decay channel at central rapidity. It is evaluated for pj > 1.3 GeV/c and |y| < 0.9, and 
extrapolated to pj > 0. The fraction of non-prompt J/ y/ in the inclusive ]/y/ yields shows no de¬ 
pendence on the charged-particle multiplicity at central rapidity. Charm and beauty hadron relative 
yields exhibit a similar increase with increasing charged-particle multiplicity. The measurements are 
compared to PYTHIA 8, EPOS 3 and percolation calculations. 
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1 Introduction 


The study of the production of hadrons containing heavy quarks, i.e. charm and beauty, in proton- 
proton (pp) collisions at the Large Hadron Collider (LHC) provides a way to test calculations based 
on perturbative Quantum Chr'omodynamics (pQCD) at the highest available collision energies. The 
inclusiveproduction cross sections of charm mesons measured in pp collisions at the LHC at both 
central lU 0] and forward |0] rapidity are described by theoretical predictions based on pQCD cal¬ 
culations with the collinear factorisation approach at next-to-leading order (e.g. in the general-mass 
variable-flavour-number scheme, GM-VLNS 0]) or at fixed order with next-to-leading-log resumma¬ 
tion (LONLL iI^-IbI]) within theoretical uncertainties. The comparisons suggest that charm production is 
under (over) estimated by the central values of the LONLL (GM-VLNS) calculations. The measured D- 
meson production cross sections in pp collisions at the LHC can also be described by pQCD calculations 
performed in the framework of kj-factorisation in the leading order (LO) approximation Jill . Beauty 
production cross section measurements in pp collisions at ^/s = 7 TeV 1 1^ 14 ] are well described by 
implementations of LONLL and GM-VLNS JiQ. In the case of B mesons, the measured cross sec¬ 
tions are close to the central value of the LONLL and GM-VLNS predictions. A similar situation was 
observed in pp collisions at ^/s = 1.96 TeV at the LNAL Tevatron collider ilhUlSlI . 


The measurement of heavy-flavour production in pp collisions as a function of the charged-paiticle mul¬ 
tiplicity produced in the collision could provide insight into the processes occurring in the collision at 
the partonic level and the interplay between the hard and soft mechanisms in particle production. These 
aspects are expected to depend on the energy and on the impact parameter (the distance between the 
colliding protons in the plane perpendicular to the beam direction) of the pp collision 1119142IH. In the 
impact parameter representation of protonproton collisions, the overlap of the nucleon wave functions in 
proton-proton collisions can be described by a geometrical picture with two separate transverse distance 
scales: the impact parameter of the collision and the transverse spatial partonic distribution 
2411. In particular, pp collisions with a hard paiton-parton scattering are predicted to be more central (i.e. 
have smaller impact parameter) than minimum-bias events i20L 12511. 


The NA27 Collaboration observed in 1988 that the average charged-particle multiplicity in events with 
open charm production was higher by about 20% than in events without charm production J3l. A 
softening of the momentum spectra of hadrons produced in association with charm was also observed. 
This result was interpreted as a consequence of the more central nature of collisions leading to charm 
production. 

At LHC energies, two additional contributions to charm production and its relation to multiplicity have 
to be considered. The first effect is the likely larger amount of gluon radiation associated to the short 
distance production processes at larger energies and _pmticle transverse momenta. The second is the con¬ 
tribution of Multiple-Parton Interactions (MPI) 127142911. i.e. several hard partonic interactions occurring 
in a single pp collision. In this context, pQCD-inspired models describe the final-state particles produced 
in hadronic collisions with a two-component approach, namely an initial hard partonic scattering pro¬ 
cess, that gives rise to collimated clusters of hadrons (jets), and an underlying event, consisting of the 
final-state particles that are not associated with the initial hard scattering. While the hard scattering pro¬ 
cess can be computed with a pQCD approach, the description of the underlying event, which is thought 
to be dominated by particles produced in soft processes and by perturbative (mini)jets with relatively 
small transverse momenta (soft MPIs), is based on a phenomenological model. In particular, pQCD- 
based models of MPIs provide a consistent way to describe high mu^licity pp collisions, and have 
been implemented in recent Monte Carlo generators like PYTHIA 6 JsCT . PYTHIA 8 |31], and HER- 
WIG 0211 . Measurements by the CMS Collaboration of jet and underlying event properties as a function 
of multiplicity in pp collisions at ^/s = 1 TeV can be better described by event generators including 
MPI 031.13411. The analysis of minijet production performed by the ALICE Collaboration 13511 indicates 
that high multiplicities in pp collisions are reached through a high number of MPIs and a higher than av- 
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erage number of fragments per parton. Upward fluctuations of the gluon density in the collidin g pro tons 
are also advocated to describe the results from high multiplicity pp collisions at the LHC 12 iL l36l. 13711 . 
Indeed, the transverse structure of the proton, as probed in hard partonic scattering processes, is predicted 
to play a crucial role in defining the underlying event structure and the probability of MPIs ll25n . In the 
heavy-flavour sector, the LHCb Collaboration reported measurements of double charm production in pp 
collisions at the LHC (D^ +X, J/\j/+X and + where A = D°, D+, D^, A+l, which suggest that 
MPIs also play a role at the hard momentum scale relevant for cc production lH, 39]. 


The ALICE Collaboration published the first measurement of inclusive J/xj/ production as a function 
of charged-particle multiplicity, expressed as the pseudo-rapidity density of charged particles dAch/dt] 
at mid-rapidity, in pp collisions at = 1 TeV 14(111 . An approximately linear increase of the yield of 
J/ i/r with the charged-particle multiplicity was observed in a multiplicity range reaching four times the 
average multiplicity (dAch/dt]). The measurements at |y| < 0.9 and 2.5 < y < 4.0 were compatible 
within the uncertainties. Both the larger amount of gluon radiation and the contribution of MPI in 
collisions where heavy quarks are produced can induce a correlation between the yield of quarkonia 
and the charged-particle multiplicity produced in the collision. The measured rise of J/i/r yield with 
increasing multiplicity can also be described in the framework of string interaction or parton saturation 
models. In particular, in Ref. 14IH a stronger-than-lineai‘ trend in the high density domain is anticipated 
as a consequence of the interaction (overlap) of strings, which reduces the effective number of sources 
for soft-particle production. The increasing trend of J/i/r yield with multiplicity is also described in a 
framework in which high multiplicities are attained in pp collisions due to the contribution of higher 
Lock states in the proton, leading to a larger number of gluons participating in the collision 13711 . 


It is also worth pointing out that the charged-paiticle densities attained in high-multiplicity pp colli¬ 
sions at the LHC are of the same order of magnitude as those measured in semi-peripheral heavy-ion 
collisions at lower centre-of-mass energies 14211 . In those heavy-ion collisions, the measured momen¬ 
tum distributions of light hadrons indicate that the system undergoes a collective expansion, which can 
be described in terms of hydrodynamics. Recent measurements in high-multiplicity p-Pb collisions at 
= 5.02 TeV ll43l - l48n and in high-multiplicity pp collisions at the LHC 114911 indicate that such a 
collective behaviour could also be at play in these systems. If charm quarks were to follow a collec¬ 
tive motion in high-multiplicity events, their momentum spectra would be altered, and the heavy-flavour 
hadron relative yields at high multiplicity would vary as a function of px 


The measurements of the px-differential prompt D meson and non-prompt J/i/r cross sections in ot col¬ 
lisions at y/s = 7 TeV with the ALICE experiment at the LHC were published in references In 

this paper, we report the measurement of the relative open heavy-flavour production yields as a function 
of the charged-particle multiplicity in pp collisions at y/s = 7 TeV. Open charm and beauty production 
is measured by reconstructing prompt D mesons and non-prompt J/i/r, i.e. J/r/r mesons coming from the 
decay of beauty hadrons. The experimental setup and the multiplicity estimation are described in Secs.|2] 
and 121 respectively. Prompt D+, D*+ mesons were measured at central rapidity, |y| < 0.5, in six 
multiplicity intervals and five pi intervals from 1 GeV/c to 20 GeV/c (Sec.lDl. The non-prompt fraction 
of J/i/r production was measured in the rapidity interval |y| < 0.9 in five multiplicity intervals and for 
px > L3 GeV/c and extrapolated to px > 0 (Sec.[5]l. The relative yields as a function of charged-particle 
multiplicity are compared in Sec.0 Linally, model calculations are discussed and compared with data in 
Sec. El 


2 Experimental apparatus and data sample 


The ALICE apparatus l5 ih consists of a central barrel detector covering the pseudo-rapidity interval 
|t]| < 0.9, a forward muon spectrometer covering the pseudo-rapidity interval —4.0 < rj < —2.5, and a 
set of detectors at forward and backward rapidities used for triggering and event characterization. In the 
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following, the subsystems that are relevant for the D meson and non-prompt J/y analyses are described. 


The central barrel detectors are located inside a large solenoidal magnet, which provides a magnetic 
field of 0.5 T along the beam direction (z axis in the ALICE reference frame). Tracking and particle 
identification are performed using the information provided by the Inner Tracking System (ITS), the Time 
Projection Chamber (TPC) and the Time Of Plight (TOP) detectors, that have full azimuthal coverage 
in the pseudo-rapidity interval |t]| < 0.9. The detector closest to the beam axis is the ITS, which is 
composed of six cylindrical layers of silicon detectors, with radial distances from the beam axis ranging 
from 3.9 cm to 43.0 cm. The two innermost layers, with average radii of 3.9 cm and 7.6 cm, are equipped 
with Silicon Pixel Detectors (SPD). The two SPD layers, covering the pseudo-rapidity ranges of |t] | < 2.0 
and IT] I < 1.4 respectively, have 1200 SPD readout chips. The two intermediate layers are made of Silicon 
Drift Detectors (SDD), while Silicon Strip Detectors (SSD) equip the two outermost layers. The high 
spatial resolution of the silicon sensors, together with the low material budget (on average 7.7% of a 
radiation length for tracks crossing the ITS perpendicularly to the detector surfaces, i.e. rj =0) and the 
small distance of the innermost layer from the beam vacuum tube, allow for the measurement of the 
track impact parameter in the transverse plane (do), i.e. the distance of closest approach of the track to 
the primary vertex in the plane transverse to the beam direction, with a resolution better than 75 /im for 
transverse momenta pj > 1 GeV/c 15211 . The SPD provides also a measurement of the multiplicity of 
charged particles produced in the collision based on track segments (tracklets) built by associating pairs 
of hits in the two SPD layers. 


At larger radii (85 < r< 247 cm), a 510 cm long cylindrical TPC 115311 provides track reconstruction with 
up to 159 three-dimensional space points per track, as well as particle identification via the measurement 
of the specific energy deposit AE/Ax in the gas. The charged particle identification capability of the TPC 
is supplemented by the TOP 15411 . which is equipped with Multi-gap Resistive Plate Chambers (MRPCs) 
located at radial distances between 377 and 399 cm from the beam axis. The overall TOP resolution 
including the uncertainty on the time at which the collision took place, and the tracking and momentum 
resolution was about 160 ps for the data-taking period considered in these analyses. 


The VO detector 15511 . used for triggering and for estimating the multiplicity of charged particles in the 
forward rapidity region, consists of two arrays of 32 scintillators each, placed around the beam vacuum 
tube on either side of the interaction region at z = —90 cm and z = -1-340 cm. The two arrays cover the 
pseudo-rapidity intervals —3.7 < rj < —1.7 and 2.8 < rj <5.1, respectively. 


The data from proton-proton (pp) collisions at a centre-of-mass energy of = 7 TeV used for the 
analyses were recorded in 2010. The data sample consists of about 314 million minimum-bias (MB) 
events, corresponding to an integrated luminosity of .ifjnt — 5 nb^^ Minimum-bias collisions were 
triggered by requiring at least one hit in either of the VO counters or in the SPD (| T] | < 2), in coincidence 
with the arrival time of proton bunches from both directions. This trigger was estimated to be sensitive 
to about 85% of the pp inelastic cross section 


To enrich the data sample with high multiplicity events, a High Multiplicity (HM) trigger based on the 
multiplicity information provided by the outer SPD layer was also used. Each readout chip of the SPD 
promptly asserts a digital pulse, called PastOR bit, on the presence of at least one firing pixel. A sample 
of about 6 million events was collected applying a selection on the minimum number of readout chips 
having asserted this digital pulse. The threshold was configured to select the 0.7% of the events with 
highest number of hits in the outer SPD layer. This HM-trigger sample (J^int — 14 nb^^) provides an 
increase of statistics by a factor of about 2.8 relative to the MB trigger for events with more than 50 
tracklets, corresponding to about six times the average multiplicity. 

Only events with interaction vertex reconstructed from tracks with a coordinate |z| < 10 cm along the 
beam line were used for the analysis. In the considered data samples, the instantaneous luminosity was 
limited to 0.6-1.2 x 10^^ cm^^s^^ by displacing the beams in the transverse plane by 3.8 times the RMS 
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of their transverse profile. In this way, the interaction probability per bunch crossing was kept in the 
range 0.04-0.08, with a probability of collision pile-up below 4% per triggered event. An algorithm to 
detect multiple interaction vertices based on SPD track segments, or tracklets, was used to further reduce 
the pile-up contribution. An event is rejected from the analysed data sample if a second interaction vertex 
is found, which has at least three associated tracklets, and is separated from the first one by more than 
0.8 cm along z. This removes about 48% of the pile-up events. The remaining pile-up contamination has 
two contributions: events with pile-up of collisions with Az < 0.8 cm and events in which the piled-up 
collisions have low-multiplicity (less than three charged particles reconstructed in the SPD). In the case 
of pile-up of collisions with small separation along z, the multiplicity estimation may be biased because 
some of the tracklets of charged particles from different interactions may be added together. According 
to simulations, the number of tracklets results to be biased when the piled-up vertices are separated 
along z by less than 0.6 cm. Combining this result with the shape of the luminous region along the beam 
direction and the maximum pile-up rate of 4%, the overall probability that two piled-up events induce a 
bias in the determination of multiplicity was found to be lower than 0.3%. The fraction of events with 
biased number of tracklets increases with increasing multiplicity and it was estimated to be below 2% at 
the highest multiplicities considered in this analysis, while the resulting bias on the measured number of 
tracklets was found to be negligible in all the multiplicity classes. 


3 Multiplicity definition and corrections 


In the present analysis, the experimental estimator of the charged-particle multiplicity is the number 
of tracklets in the interval |t]| < 1.0 (Atrackiets)- Tracklets are track segments defined by combining 
the clusters in the SPD detector with the reconstructed primary vertex position. Tracklets are required 
to point to the primary interaction vertex within ±1 cm in the transverse plane and ±3 cm in the z 
direction 0, Isih . This multiplicity estimator is the same as was used in previous studies performed for 
inclusive J/y production Idol] . Monte Carlo simulations of the detector response have shown that Atrackiets 
is proportional to the pseudo-rapidity density of the generated charged primary particles, dAch/dt], within 
2%. Primary particles are defined as prompt particles produced in the collision and all decay products, 
except products from weak decays of strange particles. The pseudo-rapidity coverage of the SPD detector 
changes with the position of the interaction vertex along the beam line, Zvtx> and with time due to the 
variation of the number of inactive channels. The detector response over the analysed data taking period 
is equalised by means of a data-based correction, which is applied on an event-by-event basis depending 
on Zvtx and time. 


The measurements in the Atrackiets £ [1^49] interval are performed using minimum-bias triggered data, 
while those in the [50,80] range exploit the SPD-based HM trigger described above. The HM trigger is 
fully efficient for events with Atrackiets > 65. The number of events and the D-meson candidate invariant 
mass distributions were corrected for the HM trigger inefficiency in the Atrackiets G [50,65] range by 
means of a data-driven re-weighting procedure. The Auackiets-dependent event weights were defined 
from the ratio of the measured distributions of the number of tracklets in the HM and minimum-bias 
trigger samples. The effect of this correction on the per-event raw yield was of about 2.5%. The average 
dAch/dt] of events in the highest Atrackiets interval was determined from the minimum-bias sample. 


The analysis results are presented as a function of the relative charged-particle multiplicity at central 
rapidity, (dAch/dq)^/(dAch/dq), where (dAch/dq) = 6.01 ±0.01 (stat.) Iqj® (syst.) is measured in in¬ 
elastic pp collisions at = 7 TeV with at least one charged particle in |t]| < 1.0 115VH . The relative 
quantities are used to minimise the experimental uncertainties and to facilitate the comparison with other 
measurements and models. The considered Atrackiets intervals and the corresponding relative charged- 
particle multiplicity ranges are summarised in Table [T] The highest Atrackiets interval considered in the 
analysis extends to a multiplicity of about 9 times the (dAch/dq) of inelastic pp collisions and the av¬ 
erage multiplicity of events in this Atrackiets interval is about six times the (dAch/dq). The uncertainty 
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on (dAch/dT])^y(dAch/dT]) is 6%; it includes the influence of (i) the determination of the Atrackiets to 
dAch/dt] proportionality factor, 2%, (ii) its possible deviation from linearity, 5%, (in) and the uncer¬ 
tainty on the measured (dAch/dt]). 


Atracklets 

(dAch/dT])l 

(dAch/dT])V(dAch/dT]) 

A°“ /lO^ 

^’events/ 

IVev/L/lO*' 

[1,8] 

2.7 

n 45+0.uj 

^•+^-0.03 

155.1 

- 

[4,8] 

3.8 


- 

89.0 

[9,13] 

7.1 

1 1 Q+0-07 
^•^^-0.07 

46.2 

50.5 

[14,19] 

10.7 

1 70+O.IO 

32.0 

35.5 

[20,30] 

15.8 

2 63+°+^ 

24.7 

28.0 

[31,49] 

24.1 

4 01 +*’■2^ 
+•'^^-0.25 

7.9 

9.5 

[50,80] 

36.7 

6 11 +°'^^ 

1.7 

- 


Table 1: Summary of the multiplicity intervals used for the analyses. The number of reconstructed track- 
lets Atrackiets, the average charged-particle multiplicity (ckVch/dTj)-', and the relative charged-particle multiplicity 
(ckVch/dTj)^/(ckVch/drj) are detailed. The number of events analysed in the various multiplicity ranges for both 
the D-meson and J/y/ analysis are reported. The number of events for the Ao-ackiets interval [50,80] are corrected 
for the high multiplicity trigger efficiency, as explained in Sec.|3] 


The analysis of D° production is also carried out as a function of the charged-particle multiplicity in 
the regions —3.7 < rj < —1.7 and 2.8 < rj < 5.1, as measured with the charge collected by the VO 
scintillator counters, Avo^ reported in units of the minimum-ionizing-particle charge. The motivation for 
studying the multiplicity dependence of charmed-meson production also with this estimator is that the 
event multiplicity and the D-meson yields are evaluated in different pseudorapidity ranges, reducing the 
effects of auto-correlations. In contrast, with the Atrackiets estimator also the D-meson decay products 
and the charged particles produced in the fragmentation of the same charm quark are included in the 
multiplicity evaluation. Monte Carlo simulations demonstrate that Avo is proportional to the charged- 
particle multiplicity in that pseudo-rapidity interval. In this paper we report D*’ relative yields as a 
function of the relative uncorrected multiplicity in the VO detector, Avo/(Avo) (see Sec. 14.4.11) . 


4 D-meson analysis 

4.1 D-meson reconstruction 


Charm production was studied by reconstructing D°, D+ and D*+ mesons, and their antiparticles, via 




their hadronic decay channels D° —)• K 71+ (with branching ratio, BR, of 3.88 ±0.05%), D+ —)• K 
(BR of 9.13 ±0.19%), and D*+ ^ D^7i+ (BR of 67.7 ±0.05%) with D° ^ K-7r+ Q. D-meson 
candidates were selected with the same strategy as described in |[ll]. The selection of D° and D+ decays 
(weak decays with mean proper decay length cT « 123 and 312 /im, respectively Js^) was based on the 
reconstruction of secondary vertices separated by few hundred microns from the interaction point. In 
the case of the D*+ strong decay, the decay topology of the produced D*’ was reconstructed. D° and D+ 
candidates were formed using pairs and triplets of tracks with the proper charge sign combination, 11] | < 
0.8, pt > 0.3 GeV/c, at least 70 associated space points (out of a maximum of 159) with ;^^/ndf < 2 of 
the momentum fit in the TPC, and at least two hits (out of 6) in the ITS, of which at least one had to be 
in either of the two SPD layers. D*+ candidates were formed by combining D° candidates with tracks 
with pj > 80 MeV/c and at least 3 hits in the ITS, out of which at least one should be in the SPD. The 
selection of tracks with | T] | <0.8 limits the D-meson acceptance in rapidity. The acceptance drops steeply 
to zero for jyj > 0.5 at low px and jyj > 0.8 at px > 5 GeV/c. A px-dependent fiducial acceptance cut. 
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Figure 1: and D+ invariant mass and D*+ mass difference distributions for selected pi and Atrackiets intervals 

for pp collisions at = 7 TeV with = 5 nb^b The D** distributions are shown for 2 < pj <A GeV /c (top- 
panels), the D+ for 4 < pj < 8 GeV /c (middle-panels), and the D*+ for 8 < px < 12 GeV/c (bottom-panels). The 
Mrackiets S [1,8], [14,19] and [31,49] intervals are shown in the left, middle and right panels respectively. The hts 
to the candidate invariant mass distributions are also shown. 


ItdI < yfid(7’T)^ was therefore applied on the D-meson rapidity. The cut value, yfid(pT)> increases from 
0.5 to 0.8 in the transverse momentum range 0 < px < 5 GeV/c according to a second-order polynomial 
function and it takes a constant value of 0.8 for px > 5 GeV/c. The selection of the decay topology was 
based on the displacement of the decay tracks from the interaction vertex, the separation between the 
secondary and primary vertices, and the pointing angle of the reconstructed D-meson momentum and its 
flight line from the primary to the secondary vertex. The selections were tuned such that a large statistical 
significance of the signal and a selection efficiency as high as possible were achieved, which resulted in 
cut values that depend on the D-meson px and species lUl]. The same selections were used in all the 
multiplicity intervals in order to minimise the effect of efficiency corrections in the ratio of the yields. 
Pion and kaon identification based on the TPC and TOE detectors were used to obtain a further reduction 
of the background. Cuts in units of resolution (at ±3 a) were applied around the expected mean values 
of energy deposit AE/Ax and time-of-flight. Tracks without TOE signal were identified using only the 
TPC information. Tracks with incompatible TPC and TOP response were considered as non-identified 
and were used in the analysis as both pion and kaon candidates. Particle identification (PID) was not 
applied to the pion tracks from the D*+ decay. This selection guarantees a reduction of the background 
by a factor of about 2 to 3 at low px, while preserving about 95% of the signal. 

The D-meson raw yields were extracted in each Atrackiets and px interval by means of a fit to the can¬ 
didate invariant mass distributions (mass difference AM = M(K7r7r) — M(K;r) for D*+). Similarly, the 
multiplicity-integrated raw yields were also evaluated for each px interval. The D*’ and D+ candidate 
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invariant mass distributions were fitted with a function composed of a Gaussian for the signal and an ex¬ 
ponential term that describes the background shape. The AM distribution of D*+ candidates, which fea¬ 
tures a narrow peak at AM ~ 145.4 MeV/c^ iSSh . was fitted with a Gaussian function for the signal and 


a threshold function multiplied by an exponential to model the background (\/AM — M;^ • 


,b{AM-M^ 




The centroids of the Gaussians were found to be compatible with the world-average masses of the D 
mesons 1581] in all multiplicity and px intervals. The widths of the Gaussian functions are indepen¬ 
dent of multiplicity and increase with increasing D-meson pj, ranging between 10 and 20 MeV/c^ for 
D° and D+ and between 600 and 900 keV/c^ for D*+ mesons, consistent with the values obtained in 
simulations. In order to reduce the influence of statistical fluctuations, the raw yields were determined 
by constraining the D-meson line shape, its mass to the world-average D-meson mass, and its width 
to the value obtained from a fit to the invariant mass distribution in the multiplicity-integrated sample, 
where the signal statistical significance is larger. Figure [U shows the D° and D+ candidate invariant 
mass distribution, and D*+ mass difference distributions, for selected pj and multiplicity intervals. The 
extraction of the raw signal yields (sum of particle and antiparticle) was possible in five pj intervals 
from 1 GeV/c to 20 GeV/c for the Mrackiets ranges reported in Table [T] The analysis covering the range 
Mrackiets G [1)49] exploited the minimum-bias triggered sample and was possible for the three D-meson 
species in three px intervals in the range between 2 and 12 GeV/c. In addition, the D° signal was 
extracted in Mrackiets G [1)49] for 1 < px < 2 GeV/c, and the D*+ signal was determined in three multi¬ 
plicity intervals for 12 < px < 20 GeV/c. The highest multiplicity interval [50,80] was studied with the 
high multiplicity triggered sample via D® mesons for 2 < px < 4 GeV/c and the three D-meson species 
for 4 < Px < 8 GeV/c. The raw yield extraction in the remaining px and multiplicity intervals for the 
different D-meson species was not performed due to the limited statistics in the analysed data sample 
and/or the large background. 


4.2 Corrections 

The yields of D mesons were evaluated for each multiplicity and px interval starting from the raw 
counts, Viaw. which were divided by the reconstruction, topological and PID selection efficiencies for 
prompt D mesons, £prompt D) and by the number of events analysed in the considered multiplicity inter¬ 
val, Advent- The results are reported as the ratio of yields in each multiplicity interval, (d^A°'*/4ydpx)^, 


to the multiplicity-integrated (average) yield, {d^N^ /dydpx). 




/dydpx 


y {d^N^° /dydp'i 


no 

raw 


\jJ pJ 
event prompt 



{K 


raw ) 


Nmb 


trigger <^tngger \<^prompt 


(^prompt D«) 


( 1 ) 


where the index j identifies the multiplicity interval. The acceptance correction, defined as the fraction of 
D mesons within a given rapidity and px interval whose decay particles are within the detector coverage, 
cancels in this ratio. D-meson raw yields have two components: the prompt D-meson contribution, and 
the feed-down contribution originating from B hadron decays. Equation[T]evaluates the yields of prompt 
D mesons under the assumption that the relative contribution to the D-meson raw yield due to the feed- 
down from B hadron decays does not depend on the multiplicity of the event, and is therefore cancelling 
in the ratio to the multiplicity-integrated values. This assumption is justified by the measurement of the 
multiplicity dependence of the B-hadron yields, via the non-prompt J/y fraction, presented in Sec.[5]and 
by PYTHIA simulations. To evaluate the yields per inelastic collisions, the number of events used for the 
normalisation of the multiplicity-integrated yield has to be corrected for the fraction of inelastic collisions 
that are not selected by the minimum-bias trigger Vmb trigger/^trigger, with Ctrigger = The 

results are also reported in Tables lATTl and IA3] without this trigger efficiency correction. It was verified 
with PYTHIA 6.4.21 lil^] Monte Carlo simulations that this minimum-bias trigger is 100% efficient for 
D mesons in the kinematic range of this measurement. 

The D-meson efficienc y co rrections were determined with Monte Carlo simulations using the PYTHIA 
6.4.21 event generator |30] with Perugia-0 tune and the GEANT3 transport code ioon . The detector 
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configuration and the LHC beam conditions were included, taking into account their evolution with time 
during the data taking period. The ^ depends on the D-meson species and on pi. For prompt D° 

mesons it is 3-4% in the 2 < px < 4 GeV/c interval and it increases up to 25-35% for px > 8 GeV/c, 
because less stringent topological selections were used at high px, where the combinatorial background 
is smaller. The efficiency for feed-down D mesons is larger by about 20-30% than for prompt D mesons. 
This is due to the fact that feed-down D mesons decay further away from the interaction vertex and are 
therefore more efficiently selected by the topological requirements. The D-meson selection efficiency 
depends also on the multiplicity of charged particles produced in the collision, because the resolution 
on the position of primary vertex improves with increasing multiplicity, providing a better resolution of 
the variables used for the topological selections. For example, the D° selection efficiency in 2 < px < 
4 GeV/c increases by about 40% from the lowest to the highest multiplicity intervals considered in this 
analysis. 


4.3 Systematic uncertainties 

Several sources of systematic uncertainty that could affect the relative yields as expressed in Eq. [T]were 
studied. Only the raw yield extraction and the feed-down subtraction contribution were found to have 
an influence on the relative yields. The influence of the raw signal extraction from the invariant mass 
distribution was evaluated by using the raw yields obtained with different approaches to separate the 
signal from the combinatorial background. The contribution to the D*’ line shape of mis-identified K 
and 71 pairs from D° decays, e.g. a D° —;r+ that passes the selection criteria as D*’ —)• ;r^K+, 
was assumed to be the same in all multiplicity intervals and was neglected in this analysis. Different 
background fit functions were considered (exponential, polynomial, linear for D° and D+; threshold, 
(AM — for D*+); the centroid and width of the Gaussians were left as free parameters in the fit 

instead of keeping them fixed to the values obtained from the multiplicity-integrated distribution; the raw 
yield was also extracted by counting the invariant mass histogram entries in a ±3a interval around the 
peak after subtracting the background evaluated by fitting the distribution side bands (i.e. excluding the 
±3a interval around the centroid). The uncertainty was estimated from the stability of the ratio of the 
raw yields £, 0 ), where the same raw yield extraction method was used in the multiplicity 

interval j and for the multiplicity-integrated result. The assigned systematic uncertainty varies from 3% 
to 15% depending on the meson species, px and multiplicity interval. 


The efficiency corrections were calculated independently for each multiplicity interval. The multiplicity 
distribution of primary charged particles in the Monte Carlo simulation, P(Ach)> was tuned to reproduce 
the measured charged-particle multiplicity 15711 . The efficiencies obtained with different Monte Carlo se¬ 
tups, that generate different initial multiplicity distributions, showed a good agreement in all multiplicity 
intervals. This effect was not considered as a source of systematic uncertainty. 


The D-meson decay tracks can be included or not in (/) the counting of the number of tracklets, resulting 
in a shift of the estimated multiplicity, and in ( ii) the determination of the primary vertex position, which 
leads to a different resolution on the vertex position and of the geometrical variables used for the D-meson 
selection. In the default configuration, the analysis was done excluding the D-meson decay tracks from 
the primary vertex determination and without excluding them from the multiplicity estimation. To check 
for possible systematic effects due to the multiplicity determination, the analysis was repeated excluding 
the D-meson decay tracks from the multiplicity estimation, obtaining compatible results. Furthermore, 
the relative yields were determined without excluding the D-meson decay tracks from the primary vertex 
determination. The influence of such variation is properly reproduced by Monte Carlo simulations, 
leading to a null effect on the corrected relative yields. Therefore this effect was not considered as a 
source of systematic uncertainty. 


The analysis was repeated for all D-meson species with different sets of topological selection criteria. 
It was verified that the corrected relative D-meson yields as defined in Eq. [T] are not sensitive to this 
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variation. This confirmed that the systematic uncertainty related to the topological selection description 
in the Monte Carlo cancels in the ratio. The influence of the PID strategy, which is based on the infor¬ 
mation of TPC and TOE detectors, was studied by also extracting the D-meson raw yields without PID 
selection criteria, which could be done only for D-meson pj > 2 GeV/ c. The ratios of the relative raw 
yields, £,o), with and without PID selections were found to be compatible with unity. As a 

consequence, this effect was not considered as a source of systematic uncertainty. 


As mentioned above, Eq.[T]describes the prompt corrected yields under the assumption that the fraction of 
prompt D mesons, /prompt> does not vary with the event multiplicity. To estimate the uncertainty related 
to this assumption, the multiplicity integrated /prompt factor was evaluated with the EONEE B-hadron 
production cross sections ||3l, the B —D-i-A decay kinematics from EvtGen |^], and the acceptance, 
reconstruction and selection efficiency of D mesons from B decays as described in |[l]]. The resulting 
/prompt values are about 85-95% depending on the D-meson pj and the applied selection criteria. The 
uncertainty due to the B feed-down contribution to the relative yields, /b = 1 — /prompt; was evaluated 
assuming a linear increase of the fraction /g/ (/b) from 1/2 to 2 from the lowest to the highest multiplic¬ 
ity interval. The resulting uncertainties vary with the pj and multiplicity range and are different for the 
three mesons. Typical values for intermediate pj at low multiplicity are and at high multiplicity 
+0 % 

-20 


4.4 Results 

The results of the D®, D+ and D*+ meson relative yields for each pj interval are presented in Eigs. |2] 
and [3] as a function of the relative charged-particle multiplicity. The relative yields are presented in the 
top panels with their statistical (vertical bars) and systematic (boxes) uncertainties except the uncertainty 
on the feed-down fraction, which is drawn separately in the bottom panels in the form of relative uncer¬ 
tainties. The position of the points on the abscissa is the average value of the relative charged-particle 
multiplicity, (dAch/dt])/ (dAch/dt]), for every Apackiets interval. The D°, D+ and D*+ meson relative 
yields are compatible in all px intervals within uncertainties. 

The average of D°, D+ and D*+ relative yields was computed for each px interval using as weights the 
inverse square of their relative statistical uncertainties. The yield extraction uncertainty was considered 
as uncorrelated systematic uncertainty. The feed-down fraction systematic uncertainty was treated as 
a correlated systematic uncertainty. The average of the D-meson relative yields for all px intervals is 
summarised in Tables lA. II and IA.21 and presented in Eig. |4(a)| The relative D-meson yields increase 
with the charged-particle multiplicity by about a factor of 15 in the range between 0.5 and six times 
(dAch/dt]). Eigure [4{b)] shows the ratios of the average of the D-meson relative yields in various px 
intervals with respect to the 2 < px < 4 GeV/c interval values. The yield enhancement is independent 
of transverse momentum within the uncertainties of the measurement. 


4.4.1 Studies with the charged-particle multiplicity at forward rapidity 

In the analysis described above, D-meson yields were measured in the same rapidity interval as the 
charged-particle multiplicity. This could lead to a bias if the particles produced in the chai'm-quark 
fragmentation and in the D-meson decay would amount to a large fraction of the measured charged 
particles. In order to study this possible bias, the measurement of the D° yields at central rapidity was 
also performed as a function of the relative charged-particle multiplicity at forward-rapidity. The charge 
collected by the VO scintillator counters, covering —3.7 < rj < —1.7 and 2.8 < rj < 5.1, was used as 
multiplicity estimator in this case. The multiplicity value Ayo was evaluated by dividing the collected 
charge by the expected average minimum-ionizing-particle charge. The D° yields were evaluated in 
intervals of Ayo; and corrected as previously described and summarised in Eq. [T] The relative yields 
of D® mesons are presented in Eig. |5] as a function of the relative mean multiplicity measured with 
the VO counters, Ayo/(Ayo). The statistical (systematic) uncertainties are represented by the vertical 
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Figure 2: D*', D+ and D*+ meson relative yields for each pj interval as a function of charged-particle multiplicity 
at central rapidity. The relative yields are presented on the top panels with their statistical (vertical bars) and 
systematic (boxes) uncertainties, except for the feed-down fraction uncertainty that is drawn separately in the 
bottom panels. D® mesons are represented by red circles, D+ by green squares, and D*+ by blue triangles. The 
position of the points on the abscissa is the average value of (dVch/dp)/ (chVch/dTj). For D+ and D*+ mesons the 
points are shifted horizontally by 1 .5% to improve the visibility. The diagonal (dashed) line is also shown to guide 
the eye. 


bars (empty boxes). Tbe systematic uncertainties due to tbe raw yield extraction and tbe B feed-down 
contribution were determined as explained in Sec. 14.31 Tbe uncertainty due to tbe unknown feed-down 
fraction evolution with tbe cbarged-particle multiplicity is drawn separately in tbe bottom panels. Tbe 
points are located on tbe x-axis at tbe average value of tbe relative mean multiplicity, Nyo / {Nyo) ■ Tbe 
uncertainty on tbe mean multiplicity values, AVo, was determined by comparing tbe mean and median 
values of tbe distributions. It was found to be below 3% for each multiplicity interval, and about 24% 
for tbe multiplicity integrated value. Tbe uncertainty on AVo/ (//vo) is not displayed on this figure. 
These results are also summarised in Tables IA.3I and IA.4I The D° relative yields increase with the 
relative uncorrected multiplicity at forward rapidity, as measured with the VO detector. The results in the 
2 < pt < 4 GeV/c and 4 < pj <S GeV/c intervals are compatible within uncertainties. The results with 
the VO multiplicity estimator indicate that the increase of the D-meson yield with the event multiplicity 
observed with the mid-rapidity estimator is not related to the fact that charmed mesons, originating 
from the fragmentation of charm quarks produced in hard partonic scattering processes, and the charged 
particle multiplicity are measured in the same pseudo-rapidity range. A qualitatively similar increasing 
trend of D-meson yield with multiplicity is indeed observed also when an rj gap is introduced between 
the regions where the D-mesons and the multiplicity are measured. 


5 Non-prompt J/ yr analysis 

5.1 Non-prompt J/ yr reconstruction 


The fraction of non-prompt J/yr in the inclusive J/yr yields, /b, was measured as a function of the 
charged-particle multiplicity by studying displaced J/yr mesons that decay into electron pairs in the ra¬ 
pidity range |y| < 0.9. This measurement, combined with the inclusive J/yr relative yield 14011 . provides 
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Figure 3: D**, D+ and D*+ meson relative yields for each px interval as a function of charged-particle multiplicity 
at central rapidity. The relative yields are presented on the top panels with their statistical (vertical bars) and 
systematic (boxes) uncertainties, except for the feed-down fraction uncertainty that is drawn separately in the 
bottom panels. D*! mesons are represented by red circles, D+ by green squares, and D*+ by blue triangles. The 
position of the points on the abscissa is the average value of (chVch/drj)/ (chVch/dTj). For D+ and D*+ mesons the 
points are shifted horizontally by 1.5% to improve the visibility. The diagonal (dashed) line is also shown to guide 
the eye. 


the multiplicity dependence of the production of beauty hadrons. J/xf/ candidates were formed by com¬ 
bining pairs of opposite-sign electron tracks. The tracks were required to have px > 1 GeV/c, at least 
70 (out of a maximum of 159) associated space points in the TPC with a ;^^/ndf of the momentum fit 
lower than 2, and to point back to the primary interaction vertex within 1 cm in the transverse plane. 
The tracks were also required to have at least one associated hit in the SPD detector, with the constraint 
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Figure 4: Average of D**, D+ and D*+ relative yields as a function of the relative charged-particle multiplicity at 
central rapidity, (a) Average of D-meson relative yields in pj intervals, (b) Ratio of the average relative yields 
in all Px intervals with respect to that of the 2 < px < 4 GeV /c interval. The results are presented in the top 
panels with their statistical (vertical bars) and systematic (boxes) uncertainties, except for the feed-down fraction 
uncertainty that is drawn separately in the bottom panels. The position of the points on the abscissa is the average 
value of (dVch/d77)/(dVch/d77). For some pj intervals the points are shifted horizontally by 1.5% to improve the 
visibility. The dashed lines are also shown to guide the eye, a diagonal on (a) and a constant on (b). 


that one of the two tracks should have a hit in the first SPD layer. Electron identification was based only 
on the TPC information. A selection of ±3a around the expected mean values of the specific energy 
deposit dE/dx for electrons was used. To further reduce the background, a ±3.5a (±3a) exclusion band 
around the expected mean specific energy deposit for pions (protons) was also applied. In order to reduce 
the combinatorial background, electron candidates compatible, together with a positron candidate, with 
being products of /-conversions (invariant mass below 100 MeV/c^) were removed. 


The measurement of /b is based on a statistical discrimination of J/ i/r mesons produced at a secondary 
vertex displaced from the primary pp collision vertex. The signed projection of the J/i/r flight distance 

onto its transverse momentum vector, pr, was constructed as L^y = (L-pr) /pt> where L is the vector 


from the primary vertex to the J/ i/r decay vertex. The pseudo-proper decay length x = (c • L^y • m) j px 
was calculated from the observed decay length using the world-average J/i/r mass m(J |'^f) = 3096.916 ± 
0.011 MeV/ iSSn . The fraction of non-prompt J/ip can be determined from a 2-dimensional un-binned 
log-likelihood fit to x and the unlike-sign di-electron invariant mass distributions. The fit procedure and 
the functions used to describe the invariant mass and the pseudo-proper decay length distributions were 
introduced in IB. 


The fraction of non-prompt J/ip as a function of the relative charged-particle multiplicity was determined 
for Px > 1.3 GeV/c in five multiplicity intervals in the Mrackiets range [4,49]. The Mrackiets G [1,3] 
range was excluded from this analysis due to the poor pseudo-proper decay length resolution, R{x), and 
the presence of a bias in the determination of x in the case of non-prompt candidates. The resolution 
of the pseudo-proper decay length is determined with Monte Caido simulations evaluating the RMS 
of the X distributions of reconstructed promptly produced J/i/r mesons. The event primary vertex can 
be computed with or without removing the decay tracks of the J/i/r candidates. The removal of the 
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Figure 5: D** meson relative yields at |y| < 0.5 for two px intervals as a function of the relative charged-particle 
multiplicity, AVo, measured at —3.7 < rj < —1.7 and 2.8 < rj < 5.1. The relative yields are presented on the 
top panels with their statistical (vertical bars) and systematic (boxes) uncertainties, except the uncertainty on the 
feed-down fraction which is drawn separately in the bottom panels. The position of the points on the abscissa is 
at the average value of AVo/ {Nyo), shifted by 1.5% to improve the visibility. The diagonal (dashed) line is also 
shown to guide the eye. 


decay tracks causes a degradation of the resolution on x, especially in the low-multiplicity intervals, as a 
consequence of the lower precision in the determination of the primary vertex with a reduced number of 
tracks. For simulated events with non-prompt J/y, the removal of the decay tracks also results in a shift 
of the primary vertex position away from the secondary decay vertex of the beauty hadrons, which is 
reflected in a systematic shift of the mean of the x distribution. However, one should consider that beauty 
quarks are always produced in pairs: the two decay tracks from the non-prompt-J/»/r, when included, 
pull the primary vertex towards the beauty hadron decay vertex, but the charged tracks from the decay 
of the second beauty quark, which enter in the barrel acceptance, pull the primary vertex in the opposite 
direction. The shift is larger in the lowest multiplicity bin where it reaches about 35 /im. This bias 
is reduced when the J/y decay tracks are kept in the evaluation of the primary vertex. The effect of 
the bias, estimated with Monte Carlo simulations, is a reduction of the measured /bQ by about 20% for 
events with Mrackiets = 4, and it becomes negligible for Atrackiets >10. Therefore, the primary vertex 
was computed considering all reconstructed tracks. To correct for the remaining bias, a modification 
in the resolution function, R{x), used to describe the non-prompt J/yr in the likelihood fit function was 
introduced, which depends on Atrackiets- In particular, the shape of the resolution function was adjusted to 
obtain a good matching between the function used to describe the non-prompt J/»/r in the likelihood fit 
(a convolution of a template of the x distribution of J/y from beauty hadron decays with the resolution 
function ll^ ) and the pseudo-proper decay length distribution of reconstructed secondary J/y from 
Monte Carlo simulations. 

Figure 0 presents the invariant mass and pseudo-proper decay length distributions for px > 1-3 GeV/c 
for each multiplicity interval together with a projection of the result of the log-likelihood fit. 

* This shift would be greater than 50 /im in the A^tracklets interval [1,3], leading to a large bias on the extracted /g value (up 
to 35%). The correction for this bias would introduce a large systematic uncertainty. 


14 



















Charm and beauty production versus charged-paiticle multiplicity 


ALICE Collaboration 



Figure 6: J/ invariant mass and pseudo-proper decay length distributions in several multiplicity intervals with 
superimposed the likelihood fit results. The contributions of the signal, the background and their sum are repre¬ 
sented with dashed, dot-dashed and full lines, respectively. In addition, the pseudo-proper decay length figures 
include the prompt and non-prompt contributions to the inclusive yields with dotted and long-dashed lines. 


5.2 Corrections 


For all multiplicity intervals, the measured fraction of non-prompt J/i/r, /g, was corrected using the 
acceptance and reconstruction efficiency of prompt, {Acc x e)prompt! and non-prompt J/i/r, (Acc x £)b, 
as 


f fi I (Accxeje \ 

® V /b (Acc x£)pi.ompty 


( 2 ) 


Here all terms refer to non-prompt J/y/ with pj > 1.3 GeV/c. The corrections for acceptance and 
efficiency were computed using Monte Carlo simulations using the GEANT3 transport code ihOn . Prompt 
J/y/ were generated with a pj distribution extrapolated from CDF measurements ilhll and a y distribution 
parameterised with the Colour Evaporation Model (CEM) 162116311 . Beauty hadrons were generated using 
the PYTHIA 6.4.21 event generator |[^ with Perugia-0 tune Ihdll . The acceptance times efficiency values 
for prompt and non-prompt J/y/ have a minimum of 8% at pT = 2 GeV/c and a broad maximum of 12% 
at pj = 7 GeV/c 16511 . The relative difference in efficiency between prompt and non-prompt J/i/r is 
only about 3%. The ratio (Acc x £)b/(Acc x £) prompt is assumed to be independent of multiplicity. The 
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uncertainty related to this assumption is discussed in the next section. 

The measured non-prompt J/i/r fractions were extrapolated from px > 1-3 GeV/c down to px = 0 using 


/r^T > 0) = a^’‘‘^-/B(pT > 1.3 GeV/c); 


^model(^T > 0) 

ymodel(^T >1.3 GcV/c) ’ 


(3) 


where represents a functional form modelled on existing data. It was calculated as the ratio of 

the differential cross section of non-prompt J/»/r, as obtained with FONLL calculations 10], to that of 
inclusive J/y, parameterised by the phenomenological function defined in 1661] : 




H2^F0NLL 


,2 _phenom 


(4) 


dydpx / dydpx 

A combined fit to the existing results of /b in pp collisions at 7 TeV E [Hi IH, in the rapidity bin 
closest to central rapidity was performed to determine the parameters of the phenomenological param- 
eterisation. The extrapolation factor obtained is = 0.991 qq]. Its uncertainties were determined by 
repeating the fit by (i) excluding the LHCb data points at forward rapidities, and (ii) using for the non¬ 
prompt J/ip cross section the upper and lower uncertainty bands of the FONLL predictions, obtained 
by varying the factorisation and renormalisation scales, instead of the central values. The uncertainties 
were determined by the maximum and minimum values obtained from these fit variations. The /b 
fractions in all multiplicity intervals were extrapolated using the same value, evaluated from the fit 
of the multiplicity integrated measurements. 


5.3 Systematic uncertainties 

The systematic uncertainty introduced by the experimental resolution on the primary vertex position was 
evaluated by repeating the fitting procedure in two alternative ways: (i) the primary vertex was evaluated 
without removing the decay tracks of the J/yr candidates. The fit was performed using the standard 
resolution function for non-prompt J/ i/r, that does not depend on multiplicity, but the x distribution of 
the non-prompt J/i/r was shifted by a multiplicity-dependent value, which was determined by the Monte 
Carlo simulation, (ii) The event primary vertex was computed after removing the decay tracks of the 
J/Y candidates and the fit was performed using the corresponding degraded resolution function R{x) and 
without any shift. The resulting uncertainties decrease with increasing multiplicity, ranging from 19% in 
the lowest multiplicity interval to 3% at the highest multiplicities. 

The uncertainty related to the extrapolation of /b from pj > 1.3 GeV/c to px > 0 was estimated with 
the method discussed above and it is about 3%. This uncertainty was assumed to be uncorrelated among 
the multiplicity intervals. 

The resolution function used in the fits is based on Monte Carlo simulations, which might introduce 
systematic effects. These were estimated by repeating the log-likelihood fits modifying the resolution 
function, R{x), according to (l/(l-|-5))./?(x/(l-|-5)), where 5 is the relative variation of the RMS of 
the resolution function, and it was varied from —0.1 to -|-0.1 to take into account the uncertainties in 
the Monte Carlo description. The systematic uncertainty due to the resolution function increases with 
multiplicity from 8% to 20%. 

The px distribution of the signal candidates (prompt and non-prompt J/y) could depend on the event 
multiplicity which could affect the shape of the resolution function which depends on the J/y Pt- The 
average px of the signal candidates was estimated from data in each multiplicity interval and found to 
be constant as a function of event multiplicity within statistical uncertainties about ±10%. The influence 
of a (px) variation on fhe resolufion funclion was defermined using Monfe Carlo simulafions: fhe px 
disfribufion was changed, considering soffer or harder px disfribufions, in order fo obfain a ±10% varia- 
fion of fhe (px)- The corresponding variafions obfained for fhe RMS of fhe resolufion function are ±7% 
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and —8.5% for the softer and harder pi distribution, respectively. The latter variations are within those 
quoted for the resolution function (±10%), therefore no additional uncertainty was included. 

The acceptance times efficiency values of prompt and non-prompt J/ i/r reconstructed for pj > 1 - 3 GeV /c 
are of the order of 10% and differ by 3%. The influence of the pj shape assumed in the simulation on 
the ratio (Acc x £)b/(Acc x £)prompt was evaluated by varying the average pt of the simulated J/i/r 
distributions within ±50%. A 1% variation in the acceptance was obtained both for prompt and non¬ 
prompt J/i/r. The corresponding variation obtained on /b through the Eq.[2]is about 1%. 


The pseudo-proper decay length shape of the combinatorial background was determined by a fit to the x 
distribution of the candidates in the sidebands of the invariant mass ilOll . By varying the fit parameters 
within their errors an envelope of distributions was obtained, whose extremes were used in the likelihood 
fit to estimate the systematic uncertainty. It increases slightly with multiplicity, ranging from 1% to 5%. 


The uncertainty on the background invariant mass shape, which was determined by fits to the invariant 
mass distributions of opposite-sign candidates in each multiplicity bin, was evaluated by using like-sign 
distributions instead, adopting the same procedure as described in lllCin . The systematic uncertainty is 
about 7%, independent of the charged-particle multiplicity. 


The shape of the x distribution of J/i/r from b-hadrons was evaluated using PYTHIA 6.4.21 1301] . The 
systematic uncertainty on its shape was computed by (i) changing the b-hadron decay kinematic, using 
EvtGen 116IH instead of PYTHA 6.4.21 or (ii) by assuming a harder and a softer b-hadron pj distribution, 
resulting in a (px) variation of about ±15%. The resulting systematic uncertainty is about 3%, constant 
with multiplicity. 


The signal invariant mass shape was fixed from fhe Monte Carlo simulation which includes the detector 
resolution effects and the radiative decays using the EvtGen 116IH package. The effect on the invariant 
mass signal shape due to the uncertainty on the detector material was studied with dedicated Monte 
Carlo simulations, where the detector material budget was varied by ±6% with respect to the nominal 
values 1^ The resulting systematic uncertainty on /b is 3% in the lowest event multiplicity interval 
and 5% in the highest one. 


The systematic uncertainties on the pseudo-proper decay length of the combinatorial background, on 
the px-exhapolation uncertainty and on the invariant mass shape of background are found or, in 
the case of assumed to be uncorrelated among multiplicity intervals. The remaining systematic 
uncertainties are (fully or partially) correlated in different multiplicity intervals. 


5.4 Results 

The relative yield of J/i/r from beauty hadron decays as a function of the charged-particle multiplicity 
was evaluated from the inclusive J/y yield and the fraction of non-prompt J/y per multiplicity interval: 


cW—"7dy dAy^/dy 
(dAj"™^"Vdy) (dfVj/v^/dy) (/b) ' 


/b is the fraction of non-prompt J/»/r in each multiplicity interval, (/b) is the fraction in the multiplicity 
integrated sample lIlOl] . and {dNy^/dy)/(dNj^^/dy) is the inclusive J/y relative yield measured for 
Pt > 0 in each multiplicity interval normalized to its value in inelastic pp collisions ll40n . All these 
quantities were measured using the same data sample and the statistical correlations were taken into 
account. In the first chai'ged-particle multiplicity class Mrackiets £ [4,8], which is used for the non-prompt 
J/Y analysis presented here, the relative yield of inclusive 3/y normalized to the inelastic cross section 
is {dN]/^/dy) j (dAj/y,/dy) = 0.41 ±0.07 (stat) ± 0.01 (syst). The values of /b extrapolated to px > 0 
were used in Eq. [5l providing the non-prompt J/ip relative yields for px > 0. The relative yields of 
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inclusive J/y were also recomputed for pj > 1.3 GeV/c and no difference was observed with respect to 
those for px > 0 within the uncertainties. 

The results for the fraction of non-prompt J/i/r for both px > 0 and px > 1.3 GeV/c, the relative yields 
of prompt and non-prompt J/ i/r in each multiplicity bin for px > 0 are summarized in Tables lA3] and lA.6l 
and shown in Fig. |7] 



Figure 7: Non-prompt J/y/ fraction as a function of the relative charged-particle multiplicity at central rapidity 
for px > 1.3 GeV/c. The vertical bars represent the statistical uncertainties, while the empty boxes stand for the 
systematic uncertainties. The width and the height of these empty boxes indicate the measurement uncertainty on 
the horizontal and vertical axis respectively. The dashed line shows the value of /b measured in the same pj range 
and integrated over multiplicity 13. The shaded area represents the statistical and systematic uncertainties on the 
multiplicity-integrated result added in quadrature. 


6 Comparison of charm and beauty production 


Figure |8(a)| presents prompt D meson and inclusive J/i/r results to compare open and hidden charm 
production. The average prompt D-meson results are shown in the 2 < px < 4 GeV/c interval with 
the px-integrated inclusive J/ ip measurement at central and forward-rapidity by the AFICE experi¬ 
ment ll40n . The results for prompt J/ip at central rapidity from this paper (px > 0) and for prompt D 
mesons (2 < px < 4 GeV/c) are compared in Fig. |8(b)| A similar increase of the relative yield with 
the charged-particle multiplicity is observed for open and hidden charm production both at central and 
forward rapidities. 


Figure |8(c)| superimposes the open charm and beauty production measurements reported in this paper 
showing the average prompt D-meson results in the 2 < px < 4 GeV/c interval and the px-integrated 
non-prompt J/i/r measurement at central rapidity. The results are compatible within the measurement 
uncertainties. 


Open charm, open beauty and hidden charm hadron relative yields present a similar increase with 

^ After the inclusive ]/¥ measurement was published in reference S, there was an improvement of the ALICE measure¬ 
ment of the inelastic cross section in pp collisions at ^/s = 1 TeV. The improved evaluation of the inelastic cross section does 
not rely on Monte Carlo, hence the systematic uncertainty is larger 1^ . To allow a proper comparison with the results reported 
here, we updated the published inclusive i/x^f measurement by the corresponding change of the trigger efficiency for inelastic 
collisions 0.864/0.85. The normalisation uncertainties were also changed from 1.5% to 


18 



























Charm and beauty production versus charged-paiticle multiplicity 


ALICE Collaboration 


Q. 

T3 

T3 


Q. 

T3 

T3 


25 

20 

15 


10 - 


0.4 

0.2 

0 

-0.2 

-0.4 


' ALICE, pp is = 7TeV 

I Average D°, D^, D*^ meson [>'|<0.5, 2<p^<4 GeV/c 
(t> J/\|/ -> eV, |y|<0.9, p^>0 
I J/V[r jiV- 2.5<y<4.0, p >0 




+6%/-3% normalization unc. not shown 
± 6% unc. on (dW/dr|) / (dW/dr|> not shown 
11111111111111111111111 

B fraction hypothesis: x 1/2 (2) at iow (high) muitipiicity 


r+4 


r+4 


(a) 


0123456789 
{dNJ6r\) / (dNJdri) 

D meson with 2 < pj <4 GeV/c and inclusive J/iff 


^.25 

Q. 

T3 

5 20 


Q. 

■Q 

■Q 


15 


10 - 


0.4 

0.2 

0 

-0.2 


' ALICE, pp i[s = 7TeV 

+ Average D^’, D’^, D*"" meson ly|<0.5, 2<p^<4 GeV/c 
+ Prompt J/y ^ eV, ly|<0.9, p >0 




+6%/-3% normalization unc. not shown 
± 6% unc. on (dW/dr|) / (dW/dr|> not shown 
11111111111111111111111 

B fraction hypothesis: x 1/2 (2) at iow (high) muitipiicity 


r+4 


r+4 


0123456789 
(dNJdri) / (dNJdri) 

(b) D meson with 2 < pj < 4 GeV/c and prompt J/t// 


^.25 

Q. 

T3 

5 20 


Q. 

■Q 

■Q 


15 


10 - 


0.4 

0.2 

0 

-0.2 


ALICE, pp is = 7Tey 

+ Average D°, D’^, D*"" meson ly|<0.5, 2<p^<4 GeV/c 
0 Non-prompt J/v eV, ly|<0.9, p >0 


* 


I 


+6%/-3% normalization unc. not shown 
± 6% unc. on (dN/dr\) / (dW/dr|> not shown 
11111111111111111111111 

B fraction hypothesis: x 1/2 (2) at iow (high) muitipiicity 


r+4 


r+4 


0123456789 
{dNJdri) / (dNJdri) 

(c) D meson with 2 < pj <4 GeV/c and non-prompt J/yr 


Figure 8: Average D meson and J/ y/ relative yields as a function of the relative charged-particle multiplicity at 
central rapidity. D-meson yields are shown for 2 < pj <4 GeV /c, while J/ y/ yields are for pj > 0. (a) Inclusive 
J/y/ results for |y| < 0.9 are represented by empty black circles lld^ . inclusive J/y/ results for 2.5 <y < 4.0 by 
black filled symbols 0, and prompt D mesons by red filled circles, (b) Prompt J/y/ results for |y| < 0.9 are 
represented by green filled crosses, and prompt D mesons by red filled circles, (c) Non-prompt J/y/ results for 
|y| < 0.9 are represented by blue filled squares, and prompt D mesons by red filled circles. The relative yields 
are presented on the top panels with their statistical (vertical bars) and systematic (boxes) uncertainties except the 
uncertainty on the feed-down fraction for D mesons, which is drawn separately on the bottom panels. The points 
are located on the x-axis at the average value of (dVch/dt])/(dVch/dtj). The diagonal (dashed) line is drawn to 
guide the eye. 


19 











































Charm and beauty production versus charged-paiticle multiplicity 


ALICE Collaboration 


charged-paiticle multiplicity. The compaitson of open and hidden heavy flavour production suggests 
that this behaviour is most likely related to the cc and bb production processes, and is not significantly 
influenced by hadronisation. The enhancement of the heavy-flavour relative yields with the charged- 
particle multiplicity is qualitatively consistent with the calculations of the contribution from MPIs to 
particle production at LHC energies |27-^. It could also be explained by the naive picture that pro¬ 
cesses with large momentum exchange might be associated to a larger amount of gluon-radiation at LHC 
energies, but no specific model implementation of this effect exists yet. The comparison of the results 
with model calculations is shown in the next section. 


7 Comparison to theoretical calculations 

Ligures |4j |7] and [8] evidence a correlation between heavy-flavour and charged-particle multiplicities. 
Heavy-flavour production is dominated by hard processes, while charged-particle yields are associated 
to the soft momentum scale processes. It is then interesting to compare our results with calculations of 
event generators, designed to be as close as possible to real events in their description of the hard and 
soft components. Even though several event gene rators are available, few of them include heavy quarks 
in a consistent way. One of these is PYTHIA ll30l .l3lh . which will be discussed in more detail in Sec. 17. II 
In Sec. 17.21 a comparison to PYTHIA 8 |31], to the EPOS 3 |71, 12] event generator results and to a 
percolation model Hlii calculation are presented. 


7.1 PYTHIA 8 simulations 

PYTHIA 8 I3II1 is the C-i-i- successor of PYTHIA 6 iBHl . One of the major improvements in PYTHIA 8 
with respect to PYTHIA 6 concerns the treatment of the MPI scenario, where the c and b quarks can 
be involved in MPI 2 —)• 2 hard subprocesses. This model improvement is fundamental for an under¬ 
standing of heavy-flavour production as a function of multiplicity, as MPI can contribute to the observed 
phenomena. Here PYTHIA 8.157 simulations with the ’’SoftQCD” process selection including colour 
reconnection and diffractive processes El are discussed, which will be referred to as PYTHIA 8. 


Heavy-flavour production in PYTHIA 8 proceeds via four main mechanisms: (i) The first (hardest) hard 
process, where the initial c/b quarks originate from the first 2 —)• 2 hard process, mostly by gluon fusion 
(gg —>• cc) or involving a c/b sea-quark (e.g. cu —)■ cu). (ii) The subsequent hard processes in MPI, 
produced via the same mechanisms as the first hard process but in consecutive interactions, that we refer 
to as hard process in MPI. Each produced gluon has a probability to split into a cc or bb pair contributing 
to heavy-flavour production. When the initial gluon originates from a hard process, either the first one or 
a subsequent process (in MPI), we refer to this process as (Hi) gluon splitting from hard process. When 
the initial gluon originates from initial or final state radiation, we refer to this process as (iv) ISR/FSR. 


The contribution of the various production processes to the total D- and B-meson production in PYTHIA 8 
for pp collisions at ^/s = 7 TeV is summarised in TablelH In the following, D mesons refer to the average 
of D°, D+, and D*+, while B mesons represent the average of B°, B+, and B*+. Initial and final state 
radiation are the main contributors to open heavy-flavour production in PYTHIA 8, corresponding to 
~ 62% for D mesons and ~ 40% for B mesons. MPI correspond to ~ 21% (~ 24%) of the D-meson 
(B-meson) production, while the first hard process is contributing ~ 11% for D mesons and ~ 36% for 
B mesons. It should be noted that in PYTHIA 8 the largest contribution to hard processes comes from c 
sea-quarks and not from gluon fusion 11741] . 


Pigure0(top panels) shows the D and B-meson production as a function of the relative charged-particle 
multiplicity calculated with PYTHIA 8. The distributions for the main production processes are shown 
independently. The top-left panel presents results for D mesons, revealing an increasing trend of the 
relative yields as a function of the relative charged-particle multiplicity for MPI, the gluon splitting 

^ In this simulation, single- and double-diffractive processes contribute to about 20% of the cross section. 
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Origin of c and b quark content 

D mesons 

B mesons 

First hard process 

11% 

36% 

gluon fusion 

2% 

15% 

c/b sea 

9% 

21% 

Hard process in MPI 

21% 

24% 

Gluon splitting from hard process 

6% 

included in ISR/FSR 

ISR/FSR 

62% 

40% 

Remnant 

< 0.2% 

< 0.4% 


Table 2: Contribution of the different production processes to the total D- and B-meson production in 
PYTHIA 8.157 IHl forpp collisions at v^= 7 TeV. 
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Figure 9: D- and B-meson relative yield as a function of the relative charged-particle multiplicity at central 
rapidity calculated with the PYTHIA 8.157 event generator M- The different c and b quark production processes 
are separated on the top panels: first hard process, hard process in multiple interactions (MPI), gluon splitting from 
hard processes and initial/final state radiation (ISR/FSR). The bottom panels present the multiplicity dependence in 
several pj intervals for prompt D-meson production, on the left for all contributions and on the right for first hard 
process only. The coloured lines represent the calculation distributions, whereas the shaded bands represent their 
statistical uncertainties at given values of (diVch/dTj) j (diVch/dTj). The diagonal (dashed) line is drawn to guide the 
eye. 


from hard processes, and the ISR/FSR contributions. This is consistent with the fact that in PYTHIA 8 
MPI and ISR/FSR contribute both to the total multiplicity and to heavy-flavour production. The first hard 
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process contribution instead shows a weaker dependence on the multiplicity: a slight increase is observed 
at low multiplicities (dAch/dt]/(dAch/dt]) < 1) followed by a saturation. The picture for B mesons, on 
the top-right panel, presents similar features as that of D mesons. The trend for the first hard process 
contribution shows an increase at low multiplicities and then saturates. The relative charged-particle 
multiplicity at which the plateau sets in is higher for B than for D mesons. The other contributions to 
particle production increase faster with multiplicity for B than for D mesons. These differences can be 
understood as being due to the larger B-meson mass, allowing a larger event activity in MPI and ISR/FSR 
processes. 

Figure |9] (bottom panels) presents the D-meson relative yields as a function of the relative charged- 
particle multiplicity in PYTHIA 8 for five pj infervals. The boffom-leff panel shows fhe frend for fhe 
sum of all confribufions, where an overall linear behaviour is observed, fhe slope of which increases wifh 
Pt- The botfom-righf panel shows fhe pj evolution for fhe firsl hai'd processes only. The relative D-meson 
yield decreases wifh mulfiplicify af low px (1 < px < 2 GeV/c), while af high px (12 < px < 20 GeV/c) 
if exhibifs a linear increase. This fealure is caused in PYTHIA 8 by fhe facf fhaf MPI are ordered by 
fheir hardness, i.e. fhe px of fhe firsf hard scattering is an upper limif for fhe subsequenf hard scafferings 
and fhe relafed ISR/FSR. Thus, charm and beaufy producfion af low px is associafed mosfly wifh low 
mulfiplicify evenfs, whereas heavy-flavour hadron producfion in high px infervals is associated fo higher 
mulfiplicify evenfs. For complefeness, fhe confribufion of MPI fo fhe fofal charged-paiticle mulfiplicify 
was sfudied. Only evenfs wifh a small number of MPI confribute fo fhe low mulfiplicify infervals, while 
high mulfiplicify evenfs are dominated by a large number of MPI, e.g. events with about five times fhe 
average mulfiplicify can have more fhan 16 parfon-parfon inferacfions in fhe even!. 

In fhe following, fhe mulfiplicify dependence for D and B-meson producfion including all confribufions 
in a given px interval, as shown in Fig.|9](bottom-leff panel) for D-mesons, is compared fo fhe measure- 
menfs. 


7.2 Comparison of data with models 


Figure [To] shows the comparison between D-meson (average of D°, D+ and D*+) production and theo¬ 
retical calculations in four px intervals. The results of the PYTHIA 8 HEl and the EPOS 3 HlTi 


event generators, and of percolation calculations 11411.17311 are represented by the red dotted line, green 
dashed or long-dashed and dotted line, and the blue dot-dashed line, respectively. The description of 
the PYTHIA 8 setup was discussed in Sec. 17.11 Figure [TT] presents px-integrated non-prompt J/»/r re¬ 
sults together with PYTHIA 8 i30l.bill calculations. The percolation model assumes that high-energy 
hadronic collisions are driven by the exchange of colour sources between the projectile and target in the 
collision Hlii. These colour sources have a finite spatial extension and can interact. In a high-density 
environment, the coherence among the sources leads to a reduction of their effective number. The source 
transverse mass determines its transverse size (oc 1/mx), and allows to distinguish between soft (light) 
and hard (heavy) sources. As a consequence, at high densities the total charged-particle multiplicity, 
which originates from soft sources, is reduced. In contrast, hard particle production is less affected due 
to the smaller transverse size of hard sources. The percolation model predicts a faster-than-linear in¬ 
crease of heavy flavour relative production with the relative charged-particle multiplicity. The D-meson 
Px-integrated percolation calculation is represented in all panels of Fig. [TO] even though in this scenario 
a Px dependence of the results is expected, such that the higher the px of the particle the stronger the 
deviation from the linear expectation. 


EPOS 3 ll7lL [720 is an event generator for various colliding systems: pp, p-A and A-A. This event 
generator imposes the same theoretical scheme in all the systems, i.e. it assumes initial conditions 
followed by a hydrodynamical evolution. Initial conditions are generated in the Gribov-Regge multiple 
scattering framework, using the ’’Parton based Gribov-Regge” formalism 17ill . Individual scatterings are 
referred to as Pomerons, and are identified with parton ladders. Each parton ladder is composed of a 
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Figure 10: Average D-meson relative yield as a function of the relative charged-particle multiplicity at central 
rapidity in different px intervals. The systematic uncertainties on the data normalisation (+6%/ — 3%), on the 
(dAch/drj)/(dAch/dp) values (±6%), and on the feed down contribution are not shown in this figure. Different 
calculations are presented; PYTHIA 8.157 ^], EPOS 3 with and without hydro 72] and a px-integrated 

percolation model . The coloured lines represent the calculation curves, whereas the shaded bands represent 

their statistical uncertainties at given values of (dAch/dp)/(dAch/dp). The diagonal (dashed) line is shown to 
guide the eye. 


pQCD hal'd process with initial and final state radiation. Non-linear effects are considered by means of a 
saturation scale. The hadronisation is performed with a string fragmentation procedure. Based on these 
initial conditions, a hydrodynamical evolution can be applied on the dense core of the collision (3-i-lD 
viscous hydrodynamics) 1721] . An evaluation within the EPOS 3 model shows that the ener gy d ensity 
reached in pp collisions at ^ = 7 TeV is high enough to apply such hydrodynamic evolution 17211 . Here 
we discuss the results of an EPOS 3.099 calculation without jet-bulk interaction, which is a process that 
produces hadrons from hard partons and quarks from the fluid. EPOS 3 without the hydro component 
(green dashed line in Pig. [TOl) predicts an approximately linear increase of D-meson production as a 
function of the charged-particle multiplicity. This linear scaling shows a px dependence, as observed in 
PYTHIA 8 with the colour reconnection scenario (red dotted line in Pig. [TOl) . although the results differ 
in magnitude. In EPOS 3, a consequence of the Parton based Gribov-Regge approach is that the number 
of MPIs is directly related to the multiplicity, i.e. Ahard process Ach Ampi- When the hydrodynamic 
evolution is considered (green long-dashed and dotted line in Pig. fTOl) . one observes a departure from a 
linear multiplicity dependence which is qualitatively comparable to that of the px-integrated percolation 
calculation. 
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Figure 11: Non-prompt J/y/ relative yield as a function of the relative charged-particle multiplicity at cen¬ 
tral rapidity for pj > 0. The systematic uncertainties on the data normalisation (-1-6%/ —3%) and on the 


(chVeh/d77)/(chVch/d77) values (±6%) are not shown in this figure. PYTHIA 8.157 13ill calculation for B mesons 
is also presented. The coloured line represents the calculation curve, whereas the shaded band represents its un¬ 
certainty at given values of (chVch/dTj)/(dAch/dTj). The diagonal (dashed) line is shown to guide the eye. 


The measurements, see Fig. [TO] and Fig. [TTJ provide evidence for an increase of the relative heavy- 
flavour yields with the relative charged-particle multiplicity which proceeds faster than linearly for high 
multiplicities. This result tends to favour calculations with a substantial deviation from linearity at high 
multiplicities such as EPOS 3.099 with hydrodynamics or the percolation model. 

8 Summary 

Charm and beauty hadron production as a function of the charged-particle multiplicity was studied in 
pp collisions at = 1 TeV. Charged-particle multiplicity at central rapidity was evaluated for events 
with at least a charged particle in the interval |t] | < 1.0. Prompt D°, D+ and D*+ meson yields were mea¬ 
sured at central rapidity (|y| < 0.5) in their hadronic decay channels in five px intervals, from 1 GeV/c 
to 20 GeV/c. The increase of the relative yield with increasing charged-particle multiplicity was found 
to be similar for D-meson species in all investigated px intervals. The average of the D°, D+ and D*+ 
relative yields increase with the relative charged-particle multiplicity faster than linearly at high multi¬ 
plicities. No px dependence is observed within the current statistical and systematic uncertainties. The 
lack of quantitative model estimates of the px dependence together with the measurement uncertain¬ 
ties prevent to conclude on a possible px dependence. A relative yield enhancement of about a factor 
of 15 with respect to the multiplicity integrated value is observed for events with six times the aver¬ 
age charged-particle multiplicity. Prompt D° relative yields were also measured as a function of the 
relative charged-particle multiplicity determined in the pseudo-rapidity intervals —3.7 < T] < —1.7 and 
2.8 < rj < 5.1. The results were found to be consistent with those obtained using the charged-particle 
multiplicity measured at central rapidity. J/y/ inclusive yields were measured earlier at central rapidity 
(|y| < 0.9) in their di-electron decay channel 0]. The non-prompt J/y/ contribution was evaluated for 
Px > 1 • 3 GeV/ c, and extrapolated to px > 0. The non-prompt J/y/ fraction does not show a dependence 
on the charged-particle multiplicity at central rapidity. 

Open charm, open beauty, and hidden charm hadron yields exhibit a similar increase with the charged- 
particle multiplicity at central rapidity. This suggests that heavy-flavour relative yields enhancement is 
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not significantly influenced by hadronisation, but more likely directly related to the cc and bb production 
processes. The heavy-flavour relative yield enhancement as a function of the charged-particle multiplicity 
is qualitatively described by: (a) PYTHIA 8.157 calculations including the MPI contributions to particle 
production percolation model estimates of the influence of colour charge exchanges during the 

interaction jdll. l73n. (c) predictions by the EPOS 3 event generator which provides a description of the 
initial conditions followed by a hydrodynamical evolution iF/ll. F^. However, the PYTHIA 8.157 13i] 
event generator seems to under-estimate the increase of heavy flavour yields with the charged-particle 
multiplicity at high multiplicities. 
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A Tables of the results 

Table IA.21 reports the results of the relative average D-meson yields per inelastic collision as a function 
of the relative charged-paiticle multiplicity at mid-rapidity in several D-meson transverse momentum 
intervals, see Fig. |4| The corresponding relative average D-meson yields normalised to the visible cross 
section instead of the inelastic one are presented in Table I A. ll 

Table IA.4I summarises the relative D** yields per inelastic collision as a function of the relative raw 
multiplicity measured with the VO detector at forward rapidity, see Fig. [5] These relative D° yields are 
also presented in Table |A3] normalised to the visible cross section. 

Table lA. 5l reports the fraction of non-prompt J/yr to the inclusive J/yr yields as a function of the relative 
charged-particle multiplicity at mid-rapidity, see Fig |7j The relative prompt and non-prompt J/»/r yields 
per inelastic collision are reported in Table IA.61 while Table |A3] presents these yields normalised to the 
visible cross section. 
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Table A.l: Average of D**, D+ and D*+ mesons relative yields for the sum of particle and antiparticle in several multiplicity and pj intervals for pp collisions at = 1 TeV 
as a function of the relative charged-particle multiplicity at central rapidity. The values are reported together with their uncertainties, which are quoted in the the order; 
statistical, systematic and feed-down contribution uncertainties. The yields reported here are not corrected by the trigger selection efficiency, they are normalised to the 
visible cross section. 
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Table A.2: Average of D**, D+ and D*+ mesons relative yields for the sum of particle and antiparticle in several multiplicity and pj intervals for pp collisions at ^/s = 1 TeV 
as a function of the relative charged-particle multiplicity at central rapidity. The values are reported together with their uncertainties, which are quoted in the the order: 
statistical, systematic and feed-down contribution uncertainties. The yields reported here are per inelastic event. 
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Table A.3: 0° meson relative yields for the sum of particle and antiparticle in several multiplicity and px intervals for pp collisions at = 7 TeV as a function of the 
relative average multiplicity in the VO detector, Vvo/ {N\q). The yields reported here are not corrected by the trigger selection efficiency, they are normalised to the visible 
cross section. 
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Table A.4: D° meson relative yields for the sum of particle and antiparticle in several multiplicity and px intervals for pp collisions at y/s = 7 TeV as a function of the 
relative average multiplicity in the VO detector, Vyo/(Avo). The yields reported here are normalised to the inelastic cross section. 
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Mracklets 

M%) 

/r(%) 

(dJv-7'/dy)/(<;77dr)x 


[4.8] 

10.1±7.8±2.5 

10.2±7.9±2.5 

0.37 ±0.07 ±0.01 

0.24±0.20^|;;JJ^ 

[9,13] 

20.8±6.9±2.7 

20.9 ±6.9 ±2.7 

0.80±0.14±0.04 

1.20 ±0.39^01^ 

[14,19] 

15.6±6.3±2.0 

15.7±6.3±2.0 

1.95 ±0.31 ±0.24 

2.06±0.84+|j;3’ 

[20,30] 

16.7 ±6.7 ±3.3 

16.8±6.7±3.3 

2.61 ±0.46 ±0.27 

2.99±1.23+gj^ 

[31,49] 

19.0±11.9±4.2 

19.0±12.0±4.2 

6.50±1.50±0.31 

8.70±5.75+|;^{ 


Table A.5: Fraction of non-prompt J/y/ measured for pj > 1 -3 GeV/c, /b(%), and extrapolated down to pj > 0, 
in the various Mrackiets intervals. Prompt and non-prompt J/y/ relative yields for pi>0 are also reported 
in the different multiplicity intervals. The first and second uncertainties correspond to the statistical and systematic 
uncertainties, respectively. The yields reported here are not corrected by the trigger selection efficiency, they are 
normalised to the visible cross section. 




0.631^;^ 

0.44 ±0.08 ±0.01 

0.28 ±0.23^^;^^ 

1 1 O+0.07 
^•^®-0.07 

0.94 ±0.17 ±0.05 

1.41±0.46;°;22 


2.29 ±0.36 ±0.28 

2.42±0.99+°;^t 

2 63^°'^^ 

^•^■^-0.17 

3.07 ±0.54 ±0.32 

3.52±1.45«;g 

^•'^^-0.25 

7.65 ±1.76 ±0.36 

10.24 ±6.76^);®® 


Table A.6: Prompt and non-prompt relative J/ y/ relative yields for pj > 0 in different multiplicity bins. The first 
and second uncertainties correspond to the statistical and systematic uncertainties respectively. The yields reported 
here are normalised to the inelastic cross section. 
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